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Abstract. Simulations of climate over the Last Millennium 
(850-1850 CE) have been incorporated into the third 
phase of the Paleoclimate Modelling Intercomparison Project 
(PM1P3). The drivers of climate over this period are chiefly 
orbital, solar, volcanic, changes in land use/land cover and 
some variation in greenhouse gas levels. While some of these 
effects can be easily defined, the reconstructions of solar, 
volcanic and land use-related forcing are more uncertain. We 
describe here the approach taken in defining the scenarios 
used in PMIP3, document the forcing reconstructions and 
discuss likely implications. 


Correspondence to: G. A. Schmidt 
(gschmidt@giss.nasa.gov) 


1 Introduction 

As part of the lead up to the Intergovernmental Panel on 
Climate Change (1PCC) Fifth Assessment report (AR5), 
modelling groups around the world are coordinating a series 
of simulations under the umbrella of the fifth Coupled 
Model Intercomparison Project (CMIP5, http://cmip-pcmdi. 
Ilnl.gov/cmip5/, Taylor et al., 2009). These simulations 
will cover the 20th century, various scenarios into the 
future and three paleo-climate simulations for the Last 
Glacial Maximum (LGM), the mid-Holocene (MH) and the 
Last Millennium (LM) in collaboration with the existing 
Paleoclimate Model Intercomparison Project Phase 111 
(PMIP3 http://pmip3.lsce.ipsl.fr/, Otto-Bliesner et al., 2009). 
This paper will discuss the choices made in determining the 
climate forcings for the LM simulations. 
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The LM simulations will start in 850 Common Era (CE) 
and go to 1850 CE where they are expected to match up and 
continue with the 1850-2005 simulations that are also being 
prescribed by CMIP5. The start date is determined by the 
lack of availability of some forcings prior to that period, and 
a desire to have the transient simulation start earlier than the 
classic Medieval period (1000-1200 CE). While the precise 
timing of the “Medieval Warm Period” or (preferably) the 
“Medieval Climate Anomaly” varies widely in the literature, 
there is strong interest in the most recent period whose mean 
temperature anomalies may have been comparable to the 
present-day (e.g., Bradley et al., 2003; Osborn and Briffa, 
2006; Mann et ah, 2009). A start at 850 CE ensures that the 
models will be able to cover all the relevant periods. 

It is likely that these experiments will be used to quantify 
the response of the models to forcings in the past in order 
to weight projections into the future (Schmidt, 2010), and 
in detection and attribution studies (Hegerl et ah, 2007). 
The experiments described below are thus intended to be 
as realistic a description of past forcings and response as 
possible, while trying hard to ensure that the uncertainties 
in the forcings are properly represented. This has lead us to 
suggest a number of alternate forcing histories for volcanic 
and solar changes rather than specifying that all models use 
a single, arbitrarily chosen, one. This approach is different 
from the one adopted for previous PMIP experiments, which 
requested that all modeling groups used similar forcing 
in order to favor model intercomparisons (Joussaume and 
Taylor, 1995; Braconnot et ah, 2007). 

The novel aspect in CM1P5 of having the same models 
and configurations used in the paleo-climate simulations 
as with the transient 20th century and future simulations 
implies that consistency - both in the overall forcings and 
in how they are imposed - with those other experiments is 
a very strong requirement. This constraint has motivated 
most of the choices described below and is one of the 
main differences between these experiments and previously 
reported simulations of this period (Jansen et ah, 2007, 
Table 6.2). 

This has a number of implications: For instance, in 

order not to bias simulations of decadal variability, solar 
reconstructions need to have a consistent representation of 
the quasi-1 1-year cycle through the record. This requires 
a synthetic cycle to be applied prior to 1610. Similarly, 
if spectral variations in solar irradiance are included in the 
20th century transient runs, then these variations need to 
be maintained over the whole period. Furthermore, many 
of the 20th century simulations will be run with interactive 
chemistry. This will allow for variations in stratospheric and 
tropospheric ozone as a function of the incoming solar UV, 
which varies both in the long term and over a solar cycle. 

Since some aspects of these reconstructions are uncon- 
strained (more so in the earlier periods), the procedure we 
adopt here is admittedly ad-hoc, though plausible. Other 
choices could be made, but these will not greatly affect 
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the range of forcings through time. Our aim is not to 
produce definitive reconstructions, but to allow the structural 
uncertainty in these reconstructions to be tested in a coherent 
environment. 

2 Forcings over the last millennium 

We describe here the first-order forcing functions that are 
relevant for this period (Table 1). This choice is determined 
by the availability of data and our current state of knowledge, 
and is therefore not anticipated to be a complete accounting. 
Other forcings, for instance related to dust or other aerosols, 
may have played a role, but cannot at this time be specified 
sufficiently for use in an intercomparison such as this. 
All timeseries and reconstructions are available via the 
Supplement and online at the PMIP3 website. 

2.1 Orbital variations 

Over the period 850 CE to present, the principal orbital 
change is the ~20day shift in the perihelion from 15 
December to 4 January, leading to a progressive increase in 
insolation in the early Northern Hemisphere summer relative 
to the later part (Berger, 1978; Laskar et al., 2004). There 
are also small decreases in eccentricity and obliquity. These 
shifts cause a decrease of about 9 W m -2 in August, and an 
increase of 3 Wm -2 in May, in high northern latitudes. In 
the Southern Hemisphere, there is a decrease of a few watts 
in November, and a similar increase in February (Fig. 1). 

2.2 Well-mixed greenhouse gases 

Variations in the principal well-mixed greenhouse gases 
(CO2, CH4 and N2O) can be derived from high resolution ice 
cores in Antarctica. Law Dome results are used exclusively 
for CO2 and CH4 (Etheridge et al., 1996, 1998; Ferretti 
et al., 2005; MacFarling Meure et al., 2006), while for 
N2O multiple cores are used in order to get a more robust 
signal (Fliickiger et al., 1999, 2002; Machida et al., 1995). 
These variations have been tied to the post-1850 record 
(Hansen and Sato, 2004) using a spline fit, after adjusting 
the CH4 record for the difference between the global mean 
and Antarctic values (16 ppb for CH4), and smoothing so that 
only decadal-scale and longer variations are retained (Joos 
and Spahni, 2008); a nominal cut-off period of 40 years was 
applied for CO2 and CH4 and of 100 years (200 years before 
1600 CE) for N 2 0 (Fig. 2). 

The variations seen are related both to natural feedbacks in 
the carbon and nitrogen cycles to changes in climate, but also 
some initial anthropogenically drive perturbations, although 
the attribution of these relatively small changes prior to the 
industrial period is difficult (Gerber et al., 2003). 
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Table 1 . Summary of Last Millennium boundary conditions. 


Feature 

PMIP3 recommendation 

Alternate solution 

Orbital parameters 

Annually varying 


Date of vernal equinox 

21 March at Noon 


Well-mixed greenhouse gases 

Annually varying 


Volcanic aerosols 

Two reconstructions (AOD, effective radius): 
GRA (Gao et al, 2008) 

CEA (Crowley et al, 2008) 


Solar irradiance 

Multiple reconstructions (including spectral variations): 
WLS (back/noback)* (Wang et al, 2005) 

1610-2005 

VSK (Vieira et al, 2010) 

850-1850 and merge to WLS (back) 

SBF (Steinhilber et al, 2009) 

850-1850 and merge to WLS (back) 

MEA (back/noback) (Muscheler et al, 2007) 
850-1610 and merge to WLS (back/noback) 

DB (back/noback) (Delaygue and Bard, 2010) 
850-1610 and merge to WLS (back/noback) 


Ozone 

Parameterisation of solar-related variations 

Same as Pl-control** 

Tropospheric aerosols 

Not prescribed 

Same as Pl-control 

Vegetation 

Land cover change (natural vegetation to crop/pasture) 

Same as Pl-control 

Ice sheets 

No changes from Pl-control 


Topography and coastlines 

Same as in Pl-control 



* The terms “back/noback” denote a solar reconstruction that does/does not include a background variation in irradiance above that driven by the 1 1 year cycle alone. 
** Pi-control is the pre-industrial CMIP5 control simulation (conditions circa 1850). 


Insolation difference (2000 - 850 CE) (W/m 2 ) 



-11 9 -7 -5 -3 -1 1 3 5 7 9 11 


Fig. 1 . Insolation difference (2000-850 CE) (Wm 2 ) as a function 
of latitude and Julian day. 


2.3 Land use/land cover (LULC) 

The reconstruction of LULC adopted by PMIP (Pongratz 
et ah, 2008) covers the time period 800 to present day. It 
is based on published maps of agricultural areas for the 
last three centuries and a country-based approach for earlier 
times that uses population data as a proxy for agricultural 
activity. The resulting dataset contains annual maps of 
the agricultural types cropland, C3 pasture, and C4 pasture 
(Fig. 3). Allocation rules were developed that determine 
which type of natural vegetation within a grid cell is affected 
by the reconstructed changes in agriculture, and extend the 
dataset to include area changes of natural vegetation (eleven 
types of forest, shrubland, and grassland). This land cover 
change has a direct impact on the albedo of the surface, but 
also indirect impacts on water cycling, surface roughness 
and soil characteristics. There is also an implied impact for 
CO 2 emissions and CO 2 uptake capacity of the terrestrial 
biosphere, but which are not considered here since they 
are implicit in the greenhouse gas concentration changes 
described above. 
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GHG Concentrations 850-1850 CE 
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Fig. 2. Well-mixed greenhouse gas concentration changes from 850 
to 1850 CE. 


Global Cropland and Pasture Area 



Fig. 3. Global mean cropland and pasture area in the Pongratz and 
Hurtt datasets. 


The reconstruction can be used through to the present 
day, with the more recent (1700 onwards) land use changes 
equivalent to those in the Ramankutty and Foley (1999) 
dataset with some regional modifications. However, for the 
simulations of the future scenarios for AR5, many modeling 
groups have adopted a recent synthesis of LULC data (Hurtt 
et al., 2009). This dataset was developed to harmonize 
historical LULC data, based on HYDE3.0 (Klein Goldewijk 
and van Drecht, 2006) from 1500 to present, with the various 
future scenarios of the Integrated Assessment Models. It 
indicates gridded area and transitions between a limited 
number of land use classes (cropland, pasture, urban area, 
primary vegetation, and secondary vegetation). The type 
of primary or secondary vegetation (forest, grassland, etc.) 


has to be determined by each climate model, e.g. from 
their dynamic vegetation modules or from maps of potential 
vegetation. Unlike for the Pongratz et al. (2008) dataset, 
no allocation rules are provided. “Secondary vegetation” 
is natural vegetation that at some point in the past has 
been affected by land use. This land use includes shifting 
cultivation and wood harvesting effects, neither of which is 
explicitly included in the Pongratz et al. (2008) dataset. Due 
to the lack of suitable input information, Hurtt et al. (2009) 
assume zero secondary vegetation in the year 1500. 

For a continuous time series of LULC data from 800 to the 
future scenarios melding the two datasets may be necessary 
for some modelling groups depending on the configuration 
and capabilities of their land model. There are natural 
discontinuities at 1500, 1700 and 1850 CE resulting from 
switches in input data sources, which might serve as good 
splice points. At the merge point, the changes in cropland 
and pasture area from the Pongratz et al. (2008) dataset can 
be scaled or applied on a grid cell-by-grid cell basis, unless 
area changes exceed the grid cell’s available agricultural or 
land area. The surplus changes in cropland and pasture 
could then be distributed on other grid cells within the same 
country. Due to the lack of wood harvest data prior to 
1700, we suggest setting the ratio of secondary vegetation to 
agricultural area constant at the merge point value for prior 
years. 

Uncertainties in the Pongratz et al. (2008) reconstruction 
can be assessed by using the maximum and minimum 
alternate reconstructions also available. 

2.4 Volcanic forcing of stratospheric aerosols 

Gao et al. (2008) (hereafter GRA) and Crowley et al. (2008) 
(hereafter CEA) provide two different reconstructions of 
aerosol optical depth due to volcanoes from stacks of well- 
dated ice core sulphate records from both polar regions. 
They differ in the transfer function for the sulphate ice 
core concentrations to aerosol optical depth (AOD) and in 
the screening processes for deciding what is and is not 
a globally important eruption. Cores from both Greenland 
and Antarctica are used to determine whether any particular 
peak is local, hemispheric or global in extent and thus 
likely to have reached the stratosphere. Historical records 
of eruptions are used to assign peaks to specific years and 
assess the likely extent of the resulting sulphates. 

2.4.1 Implementation 

The GRA procedure is to estimate the total stratospheric 
aerosol load and then use a model based estimate of the 
latitudinal, height and time distribution of the aerosols. The 
provided files contain estimates of the stratospheric aerosol 
loadings (kg SO 4 km -2 ) which can be converted into a global 
AOD by linear scaling (Stothers, 1984). The conversion to 
AOD is valid for aerosols with effective radius in the visible 
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Volcanic Stratospheric Aerosols 
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Fig. 4. The two volcanic reconstructions of aerosol optical depth (at 
550 pm) (GRA and CEA) (top), with a comparison to the modern 
estimates of Sato et al. (1993) (bottom) (note the different vertical 
scales in the two panels). 


spectral range. CEA provide an estimate of AOD directly 
from a scaling with modern observational estimates and from 
the evolution of the Pinatubo aerosols. The load can be 
estimated by reversing the above scaling. 

The difference in the two reconstructions can be seen in 
Fig. 4 (including the 20th century to allow comparison with 
the Sato et ah, 1993 data). Different events are captured in 
each case, and different amplitudes are often inferred. 

In the modern period, the dataset of Sato et ah (1993) (and 
updates) also includes variations in the effective radius (r e in 
pm) of the aerosol particles based on observations of changes 
from El Chichon (1982) and Pinatubo (1991): 

r e = 0.20 + r max (</>) x 0.75 x f(t - t Q ) 

where r max (</>) is the maximum optical depth at each latitude 
for each eruption and /(f) is a decay function (peaking at 
1.3 years and decaying to zero at 5 years). The details of this 
process may be relevant to the apparent non-linear response 
to large volcanoes in the past (Timmreck et ah, 2009), and is 
another source of uncertainty in these reconstructions. 

CEA include estimates of aerosol size, scaling to Pinatubo 
r e for eruptions up to that magnitude and for eruptions 
greater than 0.2 AOD, they use a two-thirds power scaling 
for AOD and a one-third power scaling of the AOD for r e , 
with maximum estimate of r e about 0.9 pm. For the GRA 
reconstruction, either recipe (or none) for r e may be used. 

2.5 Solar variations 

Changes in irradiance are the most important known 
component of solar activity that affects climate. Other 


aspects of solar activity - for instance, the role of 
sun’s magnetic field in shielding the planet from cosmic 
ray particles and possibly affecting clouds (Ney, 1959; 
Dickinson, 1975) are generally correlated to irradiance 
changes at least over decadal periods. 

Observations over the 30-year satellite period confirm that 
there is a quasi- 11 yr cycle in total solar irradiance (TSI) in 
phase with sunspot activity, and that the spectral signature 
of the irradiance variations is distinct, with higher energy 
frequencies (UV, X-ray etc.) varying more strongly over the 
cycle than the TSI (Lean et al., 1995; Floyd et al., 2003; 
Harder et al., 2009). Models for the irradiance have been 
successful at matching the observations based on sunspot 
numbers, magnetic activity and the area of bright faculae 
surrounding the sunspots (Frohlich and Lean, 1998; Krivova 
et al., 2003; Wenzler et al., 2006). However, there is 
some debate over whether there is a background variation 
in irradiance (the “network”) that is not easily determinable 
from observations of current activity (Foukal et al., 2006; 
Wang et al., 2005; Krivova et al., 2007), although the lower 
irradiance level during the minimum following solar cycle 
23 (i.e. in the years 2007-2009) suggests that there is indeed 
such a variation. 

Prior to the satellite era, solar activity can be tracked by 
various proxy instrumental and geophysical measurements. 
Over longer time periods, there are records of cosmogenic 
isotopes such as 14 C in tree-rings, and 10 Be or 36 C1 in 
polar ice cores. These are formed by the nuclear reactions 
between cosmic ray particles and atoms in the atmosphere 
and increase in abundance in anti-phase with solar activity. 
Connecting these proxies to the short record of satellite- 
derived solar irradiance is however difficult both technically 
and theoretically. For instance, many of the proxies do 
not overlap with the modern instruments: 14 C records have 
been contaminated by anthropogenic fossil carbon releases 
and atmospheric nuclear bomb tests, some 10 Be records 
do not have very good resolution in the 20th century and 
varying trends, and there have been measurement practice 
changes through time in sunspot counts, geomagnetic indices 
and cosmic ray intensity. The theoretical difficulties are 
also important. First, there is a small impact of historical 
geomagnetism changes on cosmic ray shielding. Second, the 
factors that affect open solar magnetic field are uncertainly 
linked to solar irradiance, though much work is being done 
to elucidate the connection (Wang et al., 2005; Krivova et al., 
2007; Vieira et al., 2010). Additionally, the relationships of 
the cosmogenic archives to global mean production levels 
may be themselves somewhat affected by climate (Field 
et al., 2009; Heikkila et al., 2008). Nonetheless, there is 
sufficient qualitative overlap of the different proxies where 
they are available to combine them quantitatively to produce 
a homogeneous record of solar activity back through time 
(Vonmoos et al., 2006; Usoskin et al., 2009). 
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2.5.1 Available input data 

The raw material out of which the solar reconstructions 
are built consists of: (i) direct satellite observations of TSI 
since 1978 (Frohlich, 2009), (ii) the Mg II index (from 
1978) (Viereck et al., 2004), (iii) Magnetograms (from 1974) 
(Wenzler et al., 2004), (iv) F10.7 radio flux (from 1947) 
(Tanaka et al., 1973), (v) instrumental records of cosmic ray 
intensity from neutron monitor counting rates since 1951 and 
ionization chambers back to 1933 (McCracken and Fleikkila, 
2003; McCracken and Beer, 2007), (vi) geomagnetic indices 
like the aa-index (from 1868) or IFIV index (from 1901) 
(Lockwood et al., 1999; Clilverd et al., 2005; Svalgaard et al., 
2004; Rouillard et al., 2007), (vii) group sunspot numbers 
(from 1610) (Hoyt and Schatten, 1998), (vii) Cosmogenic 
isotope production rates from 14 C concentrations (Muscheler 
et al., 2007; Reimer et al., 2009) (up to ~1950), and 
from multi-centennial 10 Be records from Greenland (Yiou 
et al., 1997; Beer et al., 1998; Muscheler et al., 2004; 
Berggren et al., 2009) and Antarctica (Raisbeck et al., 1990; 
Horiuchi et al., 2008). Each data source has its own 
issues with resolution, age control, measurement uncertainty, 
confounding factors and interpretation, and calibrating them 
all to the same modern measurements is complicated. 

First, the satellite measurements need to be combined into 
a consistent timeseries (Frohlich and Lean, 1998; Willson 
and Mordvinov, 2003; Dewitte et al., 2004). Unfortunately, 
there is a gap between good satellite measurements from 
mid-1989 to fall 1991 which has complicated this procedure. 
Three composite records have been published that vary 
slightly in value of the solar minimum around 1996 and 2009, 
though recent studies support the use of the Physikalisch- 
Meteorologisches Observatorium Davos (PMOD) composite 
(Krivova et al., 2009; Gray et al., 2010). Then a model 
(either statistical or physically-based) must be constructed 
for the relationship between the irradiance and other aspects 
of solar activity (Lean et al., 1995; Balmaceda et al., 2007; 
Frohlich, 2006). Using the available instrumental data and 
these models, the reconstructions can be taken back at least to 
the 19th century. Reconstructions back to 1610 are possible 
using the sunspot data using various assumption about how 
sunspots relate to the facular and network evolution. 

Wang et al. (2005) (hereafter WLS) use a physically- 
based model for the “open flux” and have produced two 
reconstructions back to 1610, one with and one without an 
independent change in the background level of irradiance 
(the “no-background” case). Alternately, Balmaceda et al. 
(2007) and Krivova et al. (2007) use a physically-based 
model of the evolution of the total solar surface magnetic 
field. Note that the WLS reconstructions (along with 
the spectral variations, Lean, 2009) are the recommended 
solar forcings for the CMIP5 20th century (1850 onwards) 
simulations. 

Before 1610, reconstructions must rely on the cosmogenic 
isotopes, principally 10 Be and 14 C. Taking past variations 
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of the geomagnetic field into account, the solar activity 
(described by the solar modulation function <J>) can be 
obtained from the isotope records. However, the relationship 
between <J> and irradiance is complicated (Usoskin et al., 
2002; Steinhilber et al., 2009). At least four groups 
have proposed <J> reconstructions from cosmogenic isotopes. 
Two of them (Muscheler et al., 2007; Delaygue and 
Bard, 2010) (hereafter MEA and DB) did not make an 
independent physical calibration to irradiance, and so we 
use a linear relation derived from the WLS modern-to- 
Maunder Minimum differences. Two other groups have 
used a chain of physically-based models to link O to 
irradiance using either the observed long-term relationship 
between total solar irradiance and open solar magnetic field 
(Frohlich, 2009; Steinhilber et al., 2009, hereafter SBF) or 
a physical modeling of the solar surface magnetic flux and its 
relationship with the isotopes (Vieira et al., 2010, hereafter 
VSK). 

2.5.2 Solar reconstructions for PMIP 

As discussed above, each reconstruction for use in a GCM 
needs to have: TSI values back to 850, a solar cycle 
(of varying amplitude) over the whole record and spectral 
variations consistent with those changes. 

We calibrate all the TSI variations to have the same value 
in the modern period. We choose to use the 1976-2006 
mean (approximately 3 cycles) from the WLS reconstruction 
(1366.14 Wm" 2 ). Recent estimates from the SORCE/TIM 
instrument suggest that the absolute value for the mean solar 
irradiance is closer to 1361 Wm -2 (Kopp et al., 2005), but 
for any models that use this or a different modern TSI (e.g., 
Mekaoui et al., 2010), we suggest simply scaling everything 
by a suitable factor. 

We define a long-term mean TSI using a 40 year triangular 
filter. A standard “11 -year” cycle is found by taking the 
maxima in the 20th century record of WLS no-background 
case, calculating the anomaly from the 40 yr smooth, 
normalising the cycle by the maximum value and averaging 
over the cycles to give a fixed 11-year shape (Fig. 5). 

Since the maximum of a solar cycle is correlated to 
the 40 year smooth TSI anomaly (Fig. 6), we can use 
this relationship to define a scaling between the long- 
term smoothed reconstruction anomaly and a plausible solar 
cycle: 

SCmax = max(0, SUmax) + 0.433 (5(f max ) - 5(f M M))/AS) 

where SC max is the solar cycle maximum anomaly at time 
fr nax, S(t) is the smoothed TSI time-series, SUmm) is the 
value of the irradiance at the Maunder Minimum (MM) 
in any specific reconstruction, and A 5 is the difference in 
smoothed TSI between the MM and modern. For values of 
TSI at or below MM values, we set the solar cycle magnitude 
to zero. Using an alternate scaling as illustrated in Fig. 6 
would increase the amplitude of the cycles by about 20%. 
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Average normalized solar cycle anomaly 



Year from Solar Max 


Fig. 5. Normalised mean solar cycle anomaly with respect to the 
40-year smoothed irradiance, derived from the no-background data 
of WLS. Greyscale is one standard deviation from the underlying 
individual cycles. 


Solar cycle maxima vs. Long term mean TS1 



Normalised 40 yr smoothed TSI anomaly 


Fig. 6. Solar cycle maxima versus the 40-year smoothed irradiance 
data from WLS. Linear regression through the origin using all local 
maxima (solid line) gives SC max =0.433 (5(f max ) — SCmmD/AS', 
(r 2 =0.47). Note that a tighter screening of which points count 
as solar cycle maxima and only using local maxima that exceed 
the long-tern mean (overlain red symbols) gives a slightly steeper 
regression of 0.53, (r 2 =0.56) (dashed line). 


Spectral solar irradiance (SSI) variations occur over solar 
cycles and the longer term, but are not purely related to the 
TSI value because of the different influences of sunspots 
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Regression of A% SSI w.r.t. A% TSI 



Fig. 7. Spectral variations of irradiance relative to the variation in 
TSI (regression of the % change in SSI at each wavelength to the % 
change in TSI). A relative magnitude of 20 implies that variations 
at that wavelength would be around 2% over a cycle, i.e. 20 x the 
~0.1% variation in TSI. The reference period used for the Lean 
(2009) and Krivova et al. (2010) curves are from the Maunder 
Minimum (1645-1715), while it is 1850 for Lean (2000). Note 
there are fewer wavelength bands in the older data. 


and faculae on different time scales at different wavelengths 
(Lean, 2000). However, linear regressions of TSI with the 
irradiance at each spectral interval explain around 97% of the 
variance at most wavelengths, implying that TSI can be used 
as a reasonable first-order predictor of the SSI prior to 1610. 
The regression of the WLS with-background TSI with the 
Lean (2009) SSI shows that the UV portions of the spectra 
have a much higher percentage variation over solar cycles 
(and presumably, longer time scales) than TSI (Fig. 7). Note 
that a number of different estimates of the spectral variations 
are emerging from new data and models (Krivova et al., 
2010; Harder et al., 2009) but are not yet available for these 
purposes. 

For the two long term reconstructions of <J>, we calibrate 
them directly to the WLS reconstructions back to 1610 
(and for both the background and no-background cases) 
and produce a consistent reconstruction from 850 to 1609 
which smoothly blends into the WLS values and the solar 
cycle phase. For the two reconstructions (VSK and SBF) 
that have independent calibrations to the Maunder Minimum 
values (one slightly greater in amplitude than the WLS with 
background case, one slightly smaller), we smoothly join 
the reconstructions in 1850 to the WLS with background 
values. We choose 1850 in these cases since the structure 
of the modern cycles is more accurate in the WLS values, 
and so that there is maximum coherence with 20th century 
transients. Other splice years are possible and could be 
examined in additional sensitivity experiments. 

Thus, we present a series of six alternative reconstructions 
(Fig. 8) that differ in their source data, processing, and 
calibration. Two of the six match up to the no-background 
WLS case, two match up with the WLS with-background 
case, and two others provide alternate estimates of the 
changes at the Maunder Minimum (matching up to WLS 
with-background in 1850). If a modelling group is using 
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Fig. 8. Multiple TSI reconstructions from 850 CE onwards, with 
the PMOD satellite-based TSI composite as reference. 


the WLS no-background case for the 20th century, they have 
a choice of two alternative continuations back to 850, while 
if they use the with-background case, they have a choice of 
four independent reconstructions for the earlier period. There 
are also four alternate estimates of TSI during the Maunder 
Minimum (and any earlier grand minima). 

The decrease in TSI at the Maunder Minimum (defined 
here as the 40-year smoothed value at 1680) compared 
to the modern average is a 0.04% decrease for the 
WLS no-background reconstructions, 0.08% for WLS with- 
background, 0.06% for SBF, and 0.10% for VSK. The 
corresponding radiative forcings are between 0.10 and 
0.23 Wm" 2 and are close to the best estimate given by IPCC 
AR4 (Forster et al., 2007, Table 2.10) of 0.19 Wm -2 . These 
changes are substantially smaller than seen in some earlier 
reconstructions (e.g., Hoyt and Schatten, 1998). However, 
the calibrations for the Maunder Minimum used previously 
are no longer justifiable in the light of new information and 
so are not considered here (Foukal et al., 2006; Forster et al., 
2007; Gray et al., 2010). 


2.5.3 Solar-related ozone changes 

In order to better match up with 20th century simulations 
that include interactive chemistry and specifically allow 
ozone to vary as a function of incoming UV, we have 
made available a parameterisation of the solar-derived ozone 
changes taken from a coupled chemistry-climate model 
(Shindell et al., 2006). This can be used in non-interactive 
runs as a modification of the basic ozone fields that is in 
phase with the solar activity and scales with the strength 
of the total solar irradiance variation relative to the modern 
mean value. 

The ozone fields come from model calculations performed 
for modern-day solar maximum and minimum conditions. 
This corresponds to a 1.1 Wm -2 irradiance change at the top 
of the atmosphere and includes spectral variations from Lean 
(2000) (Fig. 7), with the difference between equilibrium 
responses averaged over 70 years of model output. Ozone 
in the middle stratosphere increases with irradiance, with 
a maximum enhancement of ~2.5% at lOhPa. Increases are 
largest in the tropics and mid-latitudes. The magnitude of 
change decreases above and below the middle stratosphere, 
becoming negative above ~ 1 hPa. These results are similar 
to those seen in other photochemical models. Modern 
satellite observations are consistent in showing increases in 
stratospheric ozone with irradiance of 0-3%, but differ in 
both the vertical and latitudinal dependence of the increases 
(e.g., Soukharev and Hood, 2006). 

In the original modelling study (Shindell et al., 2006), 
these ozone changes increased stratospheric temperatures 
mostly around 1 hPa, with a maximum of ~ 0.6-0. 7 °C 
at mid-latitudes in both hemispheres. Warming extended 
throughout the atmospheric column, though values were 
small in the troposphere (<0.3 °C). Relative to simulations 
without ozone changes, there was an additional warming of 
0.2 °C from 100 hPa to the top of the atmosphere, and 0.05 °C 
in most of the mid-latitude and tropical troposphere. Cir- 
culation changes (associated with the strength of wintertime 
westerlies) were also directly influenced by the inclusion of 
the ozone response. The exact pattern of ozone changes 
and temperature response will however be sensitive to 
uncertainties in the spectral variations in irradiance discussed 
above. Ideally, a new parameterisation could be made for 
each substantially different spectral estimate. 


3 Discussion 

Given the variation in plausible reconstructions of the 
forcings, input data and techniques used, it is clear that 
picking a single “best” set is fraught with difficulty. Even 
though it would be easiest for modelling groups to only 
use one set of reconstructions, the very real uncertainties 
associated with that choice would not be apparent in the 
range of simulations made available for analysis. 
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Given as well the undoubted public interest in the 
solar component of post-LIA climate change, and the 
perhaps disproportionate attention given to Medieval climate 
changes, picking any one reconstruction - each of which has 
different levels of medieval-to-modern and post-LIA change 
- might be seen as attempting to bias the results of the 
models. 

We underline that these reconstructions are designed for 
a very specific purpose related to the PMIP3 simulations. 
They should not be used for sub-multi-decadal correlations to 
proxy data (since the high frequency variability is synthetic 
prior to the merge with WLS and very likely out of phase 
with the real quasi-decadal solar cycle). Neither should 
they be taken as an endorsement of the WLS, VSK or 
SBF estimates of the background changes at the Maunder 
Minimum. This calibration is still very much a focus 
of current research. Instead, we hope that the range of 
reconstructions given is wide enough so that the structural 
uncertainty in forcing can be appropriately taken into account 
in any climate change attribution studies. 

There are undoubtedly further variations that could 
be assessed. For instance, use of the Krivova et al. 
(2010) spectral changes or any new reconstruction derived 
from the Harder et al. (2009) spectral changes might 
be warranted, particularly in simulations with interactive 
ozone. Furthermore, other TSI reconstructions are also 
being developed (Shapiro et al., 2010) which suggest a larger 
change in the Maunder Minimum, but which were not 
available for use at time of writing. We hope to be able to 
provide updates to these forcings over time. 

It is worth briefly commenting on the structural impli- 
cations of the various choices being made available. It 
is relatively straightforward to calculate the envelope of 
adjusted top-of-the-atmosphere (TOA) radiative forcing for 
the various solar TSI changes (dividing by 4 and multiplying 
by the co-albedo). Forcings from greenhouse gas changes 
can be derived using formulas from Houghton et al. (2001) 
(Table 6.2) and efficacies from Hansen et al. (2005) (in order 
to account for the slight variations in physical effects of 
each change and some indirect effects such as the change 
in stratospheric water vapour as a function of changes 
in methane oxidation). We include an estimate for the 
radiative forcing associated with land use/land cover for 
surface albedo changes (Pongratz et al., 2009). This may 
be model dependent (Pitman et al., 2009) and does not 
include the implied changes in CO 2 (which are implicit in the 
greenhouse gas records), or the effects of evapotranspiration 
changes (Davin et al., 2007; Pongratz et al., 2010). Changes 
in the global mean forcing associated with orbital changes 
are small as is the net forcing due to ozone-related changes 
(Shindell et al., 2006). The AOD time series for the volcanic 
reconstructions can be used to estimate the corresponding 
radiative forcing (in units of Wm" 2 ) by multiplying it by 
—20 (Wigley et al., 2005) (though other authors suggest 
a factor of —25, Sato et al., 1993), however this does not 


Radiative Forcings (40 year smooth) (w.r.t. 850-890 CE) 
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Fig. 9. Radiative forcing from well-mixed greenhouse gases (using 
the IPCC formulas and efficacy from Hansen et al., 2005), LULC- 
related albedo change (Pongratz et al., 2009), a smoothed envelope 
of solar forcing, and the two smoothed volcanic reconstructions. 
Smoothing is with a 40-year triangular filter. 
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Fig. 10. The net radiative forcing encompassing variation from the 
solar envelope and the two volcanic reconstructions and smoothing 
as in Fig. 9. The zero-line is with reference to the earliest 40 year 
period. 


include a dependency on the r e variations, which can make 
a substantial difference for the large eruptions. 

Figure 9 shows each individual component (smoothed for 
clarity) while Fig. 10 shows the sum of all these effects 
to make a preliminary estimate of the global mean net 
effective TOA radiative forcing over time. Depending 
on the implementation of forcings and the model itself, 
the net forcing in any specific model simulation might be 
different, but this is a reasonable first order estimate of 
the likely impact on any model (c.f. Jansen et al., 2007, 
Fig. 6.13). Note that while the annual time series might be of 
use for intermediate complexity models (see Supplement), 
for GCMs we strongly suggest implementing the physical 
forcing as realistically as possible and calculating the 
resulting radiative forcing according to CMIP5 protocols. 

The most obvious feature is the difference that the volcanic 
reconstructions make on even centennial variations of the 
forcing. Given the very long term impacts of large eruptions 
(Gleckler et al., 2006), this difference is likely to be 
important in modelled responses in the 13th/14th centuries 
and in the Maunder Minimum period. Comparatively, the 
choice of solar reconstruction may not be as obviously 
important, although they have more persistent effects over 
time. Compared to the forcing histories that have been used 
previously (Jansen et al., 2007), the volcanic reconstructions 
span a similar range in radiative forcing and eruption 
frequency, while the solar forcings are slightly reduced in 
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magnitude as discussed above. Greenhouse gas and orbital 
forcings are similar to that used previously. The new forcings 
presented here (land use, solar-ozone feedbacks) have only 
rarely or never been used before in simulations of this nature. 

We anticipate that the global and hemispheric mean 
temperature changes in the model simulations will closely 
follow the smoothed forcing fields, though with amplitudes 
and lags that will vary depending on rates of ocean heat 
uptake and overall sensitivity. However, we also expect 
that the regional impact of various forcings will depend 
on the implementation of the forcing fields (whether the 
stratospheric aerosols are included as an absorbing and 
reflecting layer in the lower stratosphere, or imposed as 
a TOA radiative perturbation for instance), and the sensitivity 
of the model. While there is a higher amount of internal 
variability at the regional scale, long simulations and greater 
numbers of ensemble members will also be available for the 
analyses. Thus it is likely that the regional climate change 
fingerprints driven by the various forcings will be of most 
utility in assessing model skill (Mann et al., 2009; Schmidt, 
2010 ). 

Supplementary material related to this 
article is available online at: 
http://www.geosci-model-dev.net/4/33/2011/ 
gmd-4-33-201 l-supplement.zip. 
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